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ABSTRACT 

In this work we investigate the link between galaxy velocity dispersion, mass and other properties 
(color, morphology) with the properties of dark matter halos by comparing the clustering of galaxies 
at both fixed mass and velocity dispersion (a) . We use the Sloan Digital Sky Survey to define a volume 
limited sample of massive galaxies complete in both stellar mass (M sta r) and a with M sta r > 6 x 10 10 
h Mq and a > 75 km/s. Using this sample we show that at fixed a there is no dependence of 
the clustering amplitude on stellar or dynamical mass (Mdyn)- Conversely when M sta r or Md yn are 
fixed there is a clear dependence of the clustering amplitude on a with higher a galaxies showing a 
higher clustering amplitude. We also show that whilst when M sta r or Md yn are fixed there remains 
a dependence of clustering amplitude on morphology, there is no such dependency when a is fixed. 
However, we do see a dependence of the clustering amplitude on g — r color when both mass and a are 
fixed. Despite this, even when we restrict our samples to only ellipticals or red galaxies the relationship 
between a and clustering amplitude at fixed mass remains. It seems likely that the residual correlation 
with color is driven by satellite galaxies in massive halos being redder at fixed a. The lack of a similar 
residual morphology dependence implies that the mechanism turning satellites red is not changing 
their morphology. Our central result is that a is more closely related to the clustering amplitude 
of galaxies than either M star or Md yn - This implies that a is more tightly correlated than M star or 
Mdyn with the halo properties that determine clustering, either halo mass or age, and supports the 
notion that the star formation history of a galaxy is more closely related to its halo properties than 
its overall mass. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: halos — large-scale structure 
of universe 



1. INTRODUCTION 

In cosmological models of galaxy formation the for- 
mation history of galaxies, and hence their present day 
properties, are intimately linked with the formation his- 
tory and properties of the dark matter (DM) haloes in 
which they reside. Therefore, if we wish to understand 
the process of galaxy formation it is important to ob- 
servationally establish a link between the properties of 
galaxies and those of the halos in which they reside. 

In recent years many studies have linked galaxy prop- 
erties to halo mass using clustering, gravitational lensing, 
satellite kinematics, group catalogues or subhalo abun- 
dance matching, with the main focus being on linking 
the stellar properties (stellar luminosity, stellar mass, 
color, SFR) or those of the central black hole (lumi- 
nosity, black hole mass) to the galaxy host dark matter 
halo properties. This work has shown strong correlations 
between halo mass and stellar mass, color, star forma- 
tion rate, or morphology, such that more massive (lu- 
minous), redder (less star- forming), or more spheroidal 
galaxies live in more massive halos ( e.g. ILi et all 120061: 
Yang et all [200l [20091: iZehavi et all 120111: iMore et al l 



201H IMandelbaum et all 120061 : iLeauthaud et all 120111 ) 



For AGN the relationships are less clear although cor- 
relations do exist bet ween narrow line AGN luminos- 
ity or radio luminosity ([Wake et alll2004l : ILi et allllM 
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iWake et all l2008bl: IMandelbaum et all E009) and halo 
mass. Perhaps unsurprisingly the strongest relationship 
occurs with stellar mass, which for central galaxies in 
halos show a tight correlation with a relatively small 
scatter of 0.17 dex (jYang et all 120111) . At fixed stellar 
mass properties that indicate the star formation history 
of a galaxy then show little dependence on halo mass 
(M ore et al1l2011h . when only considering central galax- 
ies. 

In this work we focus on investigating how the dy- 
namical properties of galaxies, such as velocity disper- 
sion (a) or surface mass density (S), in addition to 
stellar mass are related to the clustering amplitude of 
galaxies and hence their host dark matter halo prop- 
erties. This is a particularly relevant question as it 
is becoming increasingly clear that many galaxy prop- 
erties such as their current star-formation rates, star- 
formation histories, metalicities, and black hole masses 
are more closely related to these dyna mical proper- 
ties than their total stellar m ass fe.g. iBender et all 
1991 iFerTarese fc MerritH [2001 iGebhardt et all 120001: 



Kauff mann et al.ll2003l : iFranx et al.ll2008l L If the forma- 



tion history and properties of the host dark matter halo 
are the major underlying cause of these galaxy proper- 
ties, we might expect that a or S are better indicators 
than stellar mass of these halo properties. 

In this paper we make use of the very large sample 
of galaxies with velocity dispersion, stellar mass, size, 
color and morphology measurements from the seventh 
data release of the Slo an Digital Sky Survey (SDSS DR7; 
lAbazaiian et al.ll2009| ). We split these galaxies into nar- 
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row bins of one parameter, in particular mass or er, and 
then see if there is any dependence of the clustering am- 
plitude on another galaxy parameter within that narrow 
bin. In this way we can investigate if there are any resid- 
ual dependencies on halo properties manifested by higher 
or lower clustering amplitudes as one parameter is varied 
whilst another remains fixed. 

Throughout this paper, we assume a flat A-dominated 
CDM cosmology with fi TO = 0.27, H a = 73km 
s _1 Mpc _1 , and a% = 0.8 unless otherwise stated. 

2. DATA 

The galaxy data used in this analysis a re gathered from 
the s eventh data release of the SDSS ()Abazaiian et al.l 
2009). We begin with the Large Scale Structure sam- 
ples of the NYU V alue Added Galaxy catalog (VAGC; 
IBlanton et al.ll2005l ). These samples have been carefully 
constructed for measurements of large scale structure and 
include corrections for fiber collisions, the tracking of 
spatially dependent completeness and appropriate ran- 
dom catalogs all of which are required to accurately mea- 
sure clustering. The sample we use has an r band mag- 
nitude range of 14.5 < r < 17.6. In addition the NYU 
VAGC gives k-correc ted (to z=0.1) absolute magnitudes 
(Blanton et al. 2003), velocity dispersion measurements 
from the Princeton Spectroscopic pipeline, and circu- 
larised sersic fits for each galaxy, all of which we make 
use of in this analysis. 

For estimates of the stellar mass we make use of the 
MPA-JHU value added catalog which provides stellar 
mass estimates bas ed on stellar population fits to the 
SPSS photometry (jKauffmann et all 120031: iSalim et all 
I2007D . The overlap between the MPA-JHU and NYU 
VAGCs is close to but to quite 100% and so we remove 
the regions where they do not match from the analy- 
sis and adjust the completeness corrections appropriately 
(see section!?]). We also remove any region of the sur- 
vey that has a spectroscopic completeness less than 70%. 
This leaves an area of 7640 deg 2 and a total sample of 
521,313 galaxies. 

The SDSS velocity dispersions are measured within 
the 3" diameter SDSS fiber. We correct to a com- 
mon aperture of one eighth of an effective radius (r e ), 
the central velocity dispersion, using t he relation op = 
P /r e ) - 066 where r ap = l."5 (ICappellari et all 
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2006) . r e is taken from the best fitting circularised sersic 
profile fit. 

In addition to investigating how clustering depends on 
stellar mass and velocity dispersion, we also include the 
dynamical mass (Md yn ) as an addition galaxy mass indi- 
cator. Estimates of galaxy stellar masses are uncertain 
due to the complex nature of the star formation history, 
stellar population synthesis mod eling, extinction law and 
initial mass function (e.g. see [Marchesini et al.l 120091 : 
iMuzzin et ail 120091: iConrov et al.ll2009bO . The dvnam- 
ical mass, which is estimated purely from the velocity 
dispersion and the size of the galaxy, will not be affected 
by these same systematic uncertainties and will likely 
have errors more similar to those of the velocity disper- 
sion. In addition there is a strong and fairly tight cor- 
relatio n between stellar m ass and dynamical mass in the 
SDSS (|Tavlor et alll2010l see also Figured]) and so both 
masses are largely tracing the same physical property of 
a galaxy. 
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Figure 1. The distribution of galaxies from the parent sample 
with M star > 6 x 1O 1O M in the M star - o plane. The colored 
regions in the top panel show the high and low a samples at 
fixed M st ar and in the bottom panel they show the high and 
low M s tar samples at fixed a. 



To estimate dynamical mass we follow iTavlor et all 
(|20Tof) where 

M dyn = K v alr e /G (1) 



with 



Ki 



73.32 



10.465 + (n-0.95) 2 



0.954 



(2) 



where n is the sersic n parameter and the gravitational 
constant G = 4.3 x 10~ 3 pc M^ 1 km 2 s -2 . 

Finally we will make use of both galaxy color and 
morphology to further parameterize our galaxy samples. 
Throughout we use g — r colors from the K-corrected 
NYU VAGC absolute magnitudes. As already stated 
these are corrected to z = 0.1; although we will refer to 
them as g — r colors they are not quite g — r at rest. The 
mor phological classificat ions we use come from Galaxy 
Zoo (jLintott et alj f2011) which provides multiple visual 
classifications for each galaxy in the SDSS spectroscopic 
sample. The parameter we use is the probability that 
a galaxy is an elliptical (P e i) corrected for redshift bias 
whereby higher redshift galaxies are more likely to be 
class ified as ellipticals du e to their smaller apparent sizes 
(see ILintott et al.1 12011. for details). 

3. SAMPLES 

We define a parent sample with stellar mass > 6 x 
10 10 M S and 0.04 < z < 0.113. The stellar mass limit and 
redshift cuts are chosen to yield the largest stellar mass 
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Figure 2. The projected correlation functions of high and low a samples at fixed Mstar (left). The projected correlation 
functions of high and low M s t ar at fixed a (right). The correlation function measurements are colored to match the selection 
regions shown in Figure [1] 



limited sample of galaxies that can be denned from the 
NYU SDSS VAGC galaxy catalog, thus maximizing the 
number of galaxies available for clustering measurements. 
This relatively high stellar mass cut is also helpful as the 
SDSS velocity dispersions are only reliable at > 75 km/s 
and thus become incomplete in a at low stellar masses. 
By inspecting the distribution of velocity dispersion in 
narrow stellar mass bins we estimate that 98% of galaxies 
with stellar mass > 6 x 10 10 M Q have velocity dispersion 
in excess of 75 km/s. Finally we remove galaxies with 
unreliable a measurements by requiring that the error in 
a be less than 10%. 

Since there is a significant scatter in the relation be- 
tween stellar mass and velocity dispersion we must make 
further cuts if we wish to define samples that are com- 
plete in both. To define the velocity dispersion com- 
pleteness limit we investigate the distribution of stellar 
mass in narrow velocity dispersion bins for the full DR7 
sample. Using these measurements we find the velocity 
dispersion that defines a complete sample of galaxies at 
the stellar mass of 6 x 10 10 AfQ. We find that for galaxies 
with a velocity dispersion of 210 km/s 85% have stellar 
masses greater than 6 x 10 10 M Q and thus we restrict our 
sample to have a grater than this limit when defining 
a limited samples. 

The aim of this work is to investigate which is more 
fundamental in determining the clustering amplitude, a 
galaxy's mass, either stellar or dynamical, or its cen- 
tral velocity dispersion. We approach this question by 
measuring the clustering where we have fixed the mass 



and allowed a to vary and visa-versa. We define a se- 
ries of samples with a narrow range in one parameter 
and then take the upper and lower quartile of the sec- 
ond parameter. We begin by defining the first velocity 
dispersion range, starting at the velocity dispersion com- 
pleteness limit (210 km/s) and cutting at a dispersion of 
242.3 km/s, the limit that contains 20% of the galaxies 
with a larger than 210 km/s. We then split that sam- 
ple by either stellar mass or dynamical mass into upper 
and lower quartiles. In a similar manner we construct 
two more samples with the same interval in a with suc- 
cessively higher a limits again splitting into highest and 
lowest quartiles in either stellar mass or dynamical mass. 
Details of these samples are given in Tables [3] and [5] in 
Appendix fA"! To make the reciprocal samples, i.e. high- 
est and lowest quartiles in a at fixed stellar mass or dy- 
namical mass, we find the appropriate starting minimum 
mass and interval that produces three samples with the 
same number of galaxies as in the a range samples. We 
then spit these into the highest and lowest quartiles of 
a. Details of these samples are given in Tables [H and 
Figure 1 shows the distribution of galaxies in the a - 
M star plane for the parent sample as well as the samples 
at fixed cr and M sta r- 

We also define a series of samples in narrow stellar mass 
ranges where we take the highest and lowest quartiles in 
stellar surface mass density (£) and their reciprocals (Ta- 
bles [5] and [7]) again ensuring that they are complete in 
both stellar mass and surface mass density. This is an 
important consistency check for our measurements. It is 



4 



Wake et al. 



232.5 < a < 921.7 
79 < a < 177.2 



11.19 < M„. 



< 11.33 



star 

253.6 < a < 565.1 
92.6 < a < 198.4 



— I — I II 
11.33 < M, 
273.7 < 
1Q8.3T < fa < 219.6 



< 11.48 



star 

< 562.7 



0.1 



10 



r (h 'Mpo) 



210 < (j < 242.3 
11.1 

10.78 < M , 

star 



< io. 9i: 



242.3 < cr < 274.6 



0.1 



3013. 



274.7 < cr < 

H- 37 < M st ar < 

10.78 < M , L < 

star 



r (h-'Mpc) 



10 



4 : 

11.911 

H.07: 



Figure 3. The ratio of the projected correlation functions between high and low velocity dispersion samples at fixed stellar mass 
(left). The ratio of the projected correlation functions between high and low stellar mass samples at fixed velocity dispersion 
(right). The best fit ratio on scales 1.6 < r p < 25.1 h -1 Mpc is shown as the dashed line. 



possible that the scatter introduced by larger errors on 
M s tar could reduce any differences in relative clustering 
amplitude compared to cr, something which should be re- 
duced by using Md yn - Whilst any errors on Md yn will be 
different and likely to be smaller than those on M s tar they 
will still be larger than the error on a due to the ad- 
ditional uncertainties in fitting the profile to determine 
R e and the sersic-n parameter. However, £, which should 
show similar trends to cr, contains both the errors on 
M st ar and on R e . Therefore, any reduction in clustering 
trends due to scatter in the measurement of the physical 
parameters should be largest for this sample. 

For each pair of samples (e.g. high and low stellar 
mass at fixed dispersion) we find that the mean of the 
fixed parameter changes very little between the samples 
split into high and low quartiles. For a fair comparison 
we compare the mass ratio to the square of the cr ratio 
since mass is proportional to cr 2 . For all samples the 
ratio of the mean of the fixed parameter varies by less 
than 5% and so should have a negligible effect on the 
clustering amplitude. Conversely the ratio of the means 
of the varying parameter covers a factor of 2 to 3.5. 

In the Appendix |A] we plot redshift, velocity, mass, 
mass density and morphology distributions for each 
sample where appropriate. As expected, since we have 
defined a volume limited parent sample and have only 
made cuts where we are complete in cr, Mg en and mass, 
the redshift distributions for all samples are the same. 



4. CLUSTERING MEASUREMENT TECHNIQUES 

The two-point correlation function, £(r), is defined as 
the excess probability above Poisson of finding an object 
at a physical separation r from another object. This is 
calculated by comparing the number of pairs as a func- 
tion of r in our galaxy catalogs with the number in a 
random catalog that covers the same volume as our data. 

Since the space density of galaxies in our subsamples 
is very low shot noise would be a significant source of er- 
ror for measurements of their auto-correlation functions. 
However, we can overcome this problem by using a much 
denser sample of galaxies to trace the underlying den- 
sity field and measure the clustering of our sub-samples 
relative to this larger sample using the cross-correlation 
function. The obvious choice for the larger sample is the 
parent sample we defined above containing all galaxies 
passing our basic selection criteria. This sample contains 
almost 65,000 galaxies and so has a space density almost 
20 times higher than our largest sub-sample. Whilst the 
individual cross-correlation functions will reflect the in- 
trinsic clustering amplitude of both the sub- and parent 
samples the ratio between any two of these cross correla- 
tions will just reflect the ratio between the clustering of 
the two sub-samples in question, with the clustering of 
the parent sample in effect being 'canceled out'. 

Measuring the line-of-sight distance using redshifts in- 
troduces distortions into the correlation function. The 
effect of these distortions can be overcome by separat- 
ing the clustering signal into contributions perpendicu- 
lar (r p ) and parallel (%) to the line-of-sight (£(r p ,Tr)). 
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Figure 4. The ratio of the projected correlation functions between high and low velocity dispersion samples at fixed dynamical 
mass (left). The ratio of the projected correlation functions between high and low dynamical mass samples at fixed velocity 
dispersion (right). The best fit ratio on scales f .6 < r p < 25.1 h~ Mpc is shown as the dashed line. 



One can then integrate over the line of sight direction to 
estimate the projected correlation function 

f°° f°° rdr£(r) 

^)= 2 { =2j^ (r2 _;,y /2 - o) 

The final expression only involves the real-space corre- 
lation function £(r) showing that w v (r v ) is not com- 
prom ised by redshift space distortions (|Davis fc Peebles! 
1983). In practice it is only possible to integrate out 
to some maximum it because £(t" p , it) is poorly known on 
very large scales resulting in additional noise being intro- 
duced to the measurement. We integrate to 60/i -1 Mpc 
which is sufficiently large to include most correlated pairs 
and gives stable results. 

We make this measurement using the lLandv fc Szalavl 
( 1993) estimator in the cross-correlation form, with 

£(7p,7t) = 

DDI ( n\ \ PR /n R \ DIR ( n R \ | ^ 

RR ynDiim/ RR \ n D J RR \nnij 

where DDI are the pair counts between the parent and 
sub-sample, DR are the pair counts between parent and 
random sample, DIR are the pair counts between sub- 
sample and random sample and RR are the pair counts 
between the random sample, all split into the r p and n 
directions. We are able to use the same random sample 
in each term as all samples are complete throughout the 
same volume, i.e. they cover the same spatial extent and 
have the same redshift distributions. We calculate the 



pair counts in a grid which is logarithmically spaced in 
r p of width 0.2 log(/i _1 Mpc) and linearly spaced in the 
7T direction of width 2/i _1 Mpc with all separations in co- 
moving coordinates. The random catalogue is based on 
the one provided by the NYU VAGC with additional cuts 
to match the exact spatial coverage of the MPA/JHU 
stellar masses and redshifts matching the redshift distri- 
bution of our samples. The resulting random catalogue 
contains almost 1.9 million randoms, close to 29 times as 
many as in our parent galaxy catalog, which is a suffi- 
cient number such that our errors are never dominated 
by shot noise in the random pair counts. 

Since the individual bins in a correlation function mea- 
surement are highly correlated we need to generate accu- 
rate covariance matrices if we wish to compare our clus- 
tering measurements in a meaningful way. The simplest 
approach and the one we chose to follo w is to use jack- 
knife resampling (jScranton et al.l l2002T l . There is some 
debate about the accuracy of co yariance matric e s gen - 
erated in this way; for instance iNorberg et al.l (2009) 
suggest that the structure of jackknife covariance ma- 
trices may not be entirely accurate for w p (r p ) mea- 
surements, a l though the varian ce is reliable. Conversely 
iZehavi et afl (|2002l 12004 l2005f ) find that they can closely 
reproduce the structure and amplitude of covariance ma- 
trices generated by mock catalogs using jackknife resam- 
pling. Because of this uncertainty, it might be prefer- 
able to generate covariance matrices from large numbers 
of mock galaxy catalogues generated by populating dark 
matter halo catalogs so as to reproduce the clustering and 
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Figure 5. The ratio of the projected correlation functions between high and low surface mass density samples at fixed stellar 
mass (left) and high and low stellar mass samples at fixed surface mass density (right). The best fit ratio on scales f .6 < r p < 25.1 
h _1 Mpc is shown as the dashed line. 



density of each sample. However, since we have so many 
sub-samples, each with a complex selection, it would be 
enormously challenging to undertake such a scheme. We 
have therefore chosen to use jackknife resampling which 
should be sufficient for our purposes. We split the SDSS 
area into 146 equal area regions and then repeatedly cal- 
culate w p (r p ) removing one area at a time. These 146 
w p (r p ) measurements are then used to generate a full 
covariance matrix. 

5. RESULTS 

5.1. Is galaxy mass or a more closely related to 
clustering amplitude? 

Figure [2] shows the w p (r p ) measurements for the high 
and low a samples at fixed M sta r (left) and the high and 
low M s tar samples at fixed a (right). We find that when 
the mass is held fixed the high a samples have a higher 
clustering amplitude than the low er samples. This is true 
on all scales and for all mass ranges. However, when a is 
held fixed there appears to be little or no dependence of 



the clustering amplitude on M sta r- This is the central 
result of this paper: clustering amplitude depends on 
a at fixed mass, but does not depend on mass at fixed a. 

This result is illustrated more clearly in Figure [3] which 
shows the ratio of w p (r p ) between the high and low 
cr samples at fixed M sta r (left) and the high and low 
M s tar samples at fixed a (right). Figure 0] is similar to 
Figure [3j except M s tar is replaced by Md yn and shows 
exactly the same trends as with M s tar ■ We also show 
in Figure [5] the ratio of w p (r p ) between high and low 
E samples at fixed M sta r and high and low M sta r samples 
at fixed E. This is, in effect, the same as splitting by 
size at fixed mass or mass at fixed size. Since Md yn and 
Mgtar are highly correlated one would expect galaxies 
with smaller R e at a fixed M s tar, thus higher S, to have 
a higher a. These samples should then be analogous to 
those with varying a at fixed M s tar and indeed they show 
the same trend of higher clustering amplitude at higher 
£ when M s tar is fixed and very little dependence with 
Mstar when E is fixed. As discussed in Section [3] this also 
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implies that these observed trends are not the result of 
larger measurement errors on M sta r or Md yn compared to 
a, which could have reduced the clustering dependence 
on mass compared to a . 

To quantify these ratios in w p (r p ), which constitutes a 
ratio in galaxy bias, we find the best fitting scale inde- 
pendent amplitude on scales 1.6 < r p < 25.1 h _1 Mpc. 
These scales are sufficiently large as to be dominated 
by pair counts between DM halos, rather than galaxies 
within halos, and so should show a scale independence 



with respect to the dark matter clustering and give an 
indication of the relative linear bias of the two samples. 
For the w p (r p ) ratios shown in Figures |4] and [5] this 
appears to be the case and indeed a scale independent 
ratio is an acceptable fit for all but one of the samples. 
These fits are made using the full covariance matrices 
for the w p (rp) ratio, and are shown as the dashed lines 
in Figures [3J [4] and [5j Details of the fits are given in 
Table Q] and the best fit bias ratios and errors are plot- 
ted for all these samples in Figure |U In addition to the 
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Table 1 

Fits to correlation funtion ratios 

Sample Best fit ratio Best fit \ A Ratio = 1 x 2 Ratio = 1 prob 
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Figure 7. The relationship between dynamical and stellar 
mass for all DR7 galaxies with 0.04 < z < 0.113 and a > 
75 km/s (black contours) and those galaxies classified as ei- 
ther ellipticals (red contours) or spirals (blue contours). The 
points show the median M s tar in bins of Md yn and the solid 
lines the biweight best fit. The dotted line shows the one-to- 
one relation. This relationship is almost identical indepen- 
dent of morphology over the whole mass range covered by 
this plot, and becomes indistinguishable for the mass range 
of interested in this paper. 



best fit we also calculate the \ 2 for a w p (r p ) ratio of one 
over the same scales, and give the \ 2 value as well as the 
probability of an acceptable fit in Table [T] In all cases 
the best fit ratio is higher when a or S is varied at fixed 
mass than when a or E are fixed and the mass allowed to 
vary. Similarly when a or E are fixed the uu p (r p ) ratios 
between the high and low mass samples are consistent 
with one in all but one case whereas the varying a or 
E samples at fixed mass are nearly all inconsistent with 
a w p (r p ) ratio of one. For each parameter pair we can 
combine the samples and determine the probability that 
all are consistent with a w p (r p ) ratio of 1. Again we find 
that at fixed a or E the samples are indeed consistent 
with showing no dependence of the clustering amplitude 
on mass at the 66%, 87% and 66% levels, whereas there 
is only a 0.9%, 0.03% and 0.7% chance that all three of 
the varying a or E samples at fixed mass are consistent 
with no variation in clustering amplitude on large scales. 

On small scales the differences are just as clear with 
the varying a or E samples typically showing even larger 
w p (r p ) ratios than at large scales. The samples with 
varying mass at fixed a or E occasionally show moder- 
ately larger w p (r p ) ratios on small scales than on large 
scales, but often there is no difference with the w p (r p ) ra- 
tio remaining consistent with one on all scales. 

Overall these measurements make it clear that veloc- 
ity dispersion or stellar mass surface density are more 
closely related to galaxy bias than either stellar mass or 
dynamical mass. 
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5.2. The Importance of Morphology and Color 

We have included the distribution in elliptical proba- 
bility for each of the samples in the plots in Appendix 
[Al as it demonstrates one area that is both potentially a 
concern but also scientifically important for this analy- 
sis, a is much better correlated with morphology or color 
than stellar mass is. £ is also a better discriminator than 
mass, but not as good as a (see Wake et al. 2012 for a 
detailed study of these trends). Because of this when 
we fix a we come pretty close to fixing morphology re- 
gardless of the mass and when we split into high and low 
a samples at fixed mass we are pretty close to splitting 
into disks and spheriods (we also see a splitting by color 
in the same way). This is of course the heart of our in- 
vestigation: is a tighter correlation between a and halo 
properties (compared to the correlation between mass 
and halo properties) the explanation as to why a is a 
better indicator than stellar mass of a galaxy's stellar 
population? This does however raise a potential worry 
regarding systematics: are we measuring the same thing 
when we measure a for a disk galaxy as compared to a 
spheroidal galaxy? 

We tackle this question in several ways. Firstly we 
show in Figure [7] the observed relationship between 
M sta r and Md yn for all the galaxies in SDSS with 0.04 < 
z < 0.113 and a > 75 km/s. We further split this sample 
into spheriods or disks based on the Galaxy Zoo classifi- 
cations. The relationship is essentially identical regard- 
less of morphology over practically the entire mass range. 
There is a small deviation at lower masses but this is be- 
low the mass range considered in this work. This gives 
us confidence that the velocity dispersion is measuring 
the same physical quantity for both disks and ellipticals 
within the SDSS fiber radius. 

The second approach is a more direct one: we can see 
if morphology itself has any effect on the clustering am- 
plitude of our galaxies at fixed mass or a. We thus define 
samples with high and low elliptical probabilities in nar- 
row ranges of Md yn and a. We note that we are unable 
to just take the quartiles of the distribution in elliptical 
probability as we have done for parameters previously 
as there is a redshift dependence to the probabilities as- 
signed. This results from galaxies being harder to classify 
at higher redshift due to their smaller angular size and 
surface brightness dimming. Whilst there has been some 
correction for this effect applied to the Galaxy Zoo prob- 
abilities we find that it does not entirely remove this bias 
resulting in the samples split into upper and lower quar- 
tiles in P e i having different redshift distributions. We 
thus define our samples making use of the expectation 
that the morphological mix does not vary significantly 
over our narrow redshift range. We find the best fit lin- 
ear relations between elliptical probability and redshift 
that produce the 25% most and least likely ellipticals at 
each redshift and use these relations to define our sam- 
ples. Details of these samples are given in Tables [9] and 
1101 and their distributions are plotted in Appendix lAl 

Figure [8] shows the large scale w p (r p ) ratios of these 
high and low P i samples at fixed mass and a. As before 
we fit the w p (r p ) ratio (shown in Figure [T7] in Appendix 
[A"|l between 1.6 < r p < 25.1 h~ 1 Mpc and determine the 
probability that the ratio is consistent with one. Details 
of these fits are give in Table [2] At fixed a there is no 



evidence for any dependence of the clustering amplitude 
on morphology at any scale. At fixed mass the ratios are 
all non-zero with galaxies with a higher P c i showing a 
higher clustering amplitude than those with a low P e i- 
However, the significance of these positive w p (r p ) ratios 
is low with all being consistent with no difference in the 
clustering amplitude. This implies that the correlation of 
morphology with dispersion is not the determining factor 
in the strong correlation of clustering amplitude with a at 
fixed mass. 

Whilst the strong correlation of color with a is unlikely 
to introduce any systematics in our measurements it is 
still informative to investigate whether there is any resid- 
ual clustering dependence on color at either fixed mass 
or a . Since color is better correlated with a than mass 
one might imagine that any residual clustering depen- 
dence on color would be lower at fixed a than at fixed 
mass if the halo properties are key to determining the 
color, much as is seen with morphology above. To test 
this we show in Figure [8] the best fit large scale correla- 
tion function ratios of high and low g — r color samples 
at fixed mass and dispersion and list the fitting results 
in Table [2] (the full w p (r p ) ratios are shown in Figure 
[T8"lin Appendix [A]). Since the color evolves with redshift 
we have again had to split the samples in a similar way 
to the elliptical probability samples such that we have 
the reddest and bluest quartiles at any redshift. Now, 
unlike with the elliptical probability, we see a residual 
clustering dependence on color at both fixed mass and 
a such that red galaxies cluster more strongly than blue. 
This is most likely the result of the truncation of star 
formation in satellite galaxies in massive dark halos, re- 
sulting in satellite galaxies at fixed mass or a having 
redder c olors than centra l galaxies with the same mass 
or a (e.g. lYang et al.ll2009[ ). The fact that we see a signif- 
icant residual color dependence and not much of a resid- 
ual morphology dependence suggests that the environ- 
mental mechanism for transforming the color of a galaxy 
from red to blue does not simultaneously ch ange the mor- 
pholog y. A similar result was observed by iSkibba et al.1 
(2009) using the marked correlation function applied to 
SDSS galaxies. 

Since we do see some residual color dependence and 
perhaps a hint of morphology dependence in the cluster- 
ing amplitude at fixed mass it is worth asking if the a de- 
pendence remains when morphology or color is fixed. We 
show in Figure [9] the large scale w p (r p ) ratio of high and 
low a samples at fixed Md yn and high and low Md yn sam- 
ples at fixed a where we have restricted the galaxies to 
be either red (g — r > 0.9) or have high elliptical prob- 
abilities (P e i > 0.6). The fit details given in Table [2] 
and the full uv,(r p ) ratios shown in Figures [20] and [19] 
in Appendix [A] Despite the sample sizes being reduced 
and hence the errors increasing, the original trend of high 
a galaxies having a higher clustering amplitude at fixed 
mass remains significant and the large scale ratios stay es- 
sentially the same within the errors. There is a hint that 
the w p (r p ) ratios are slightly reduced, but by removing 
the blue or spiral galaxies we have removed many lower 
a galaxies and reduced the difference between the high 
and low a samples. This is a clear demonstration that 
galaxies with higher a cluster more strongly than galax- 
ies with lower a regardless of their mass, morphology or 
color. 
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Table 2 

Fits to correlation funtion ratios for sample split by morphology and color 



Sample 


Best fit ratio 


rSest nt \ 


rtatio = 1 X 


Ratio = 1 prob 


P c l at cr 1 


0.92 
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0.05 


2.44 


4.55 
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Pel at cr 2 


1.06 
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0.09 
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0.63 
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6.16 
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1.20 
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2.02 


0.847 


T» i Ti ,r All 

Pel at M dyn All 










13.95 


0.53 


g — r at a 1 


1.35 
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0.11 


3.35 


11.81 


0.037 


g — r at cr 2 


1.36 
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g - r at M dyn 1 


1.34 
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6.66 


16.34 
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1.50 
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12.76 
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0.11 
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34.29 
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1.29 
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13.46 
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1.48 


± 
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6. DISCUSSION 

We have shown that at fixed mass (stellar or dynam- 
ical) galaxies with high a (or S) cluster more strongly 
than galaxies with low a . Conversely, at fixed a (or S) 
our measurements are consistent with high and low mass 
galaxies clustering the same. If the clustering of galax- 
ies is dominated by the clustering of the DM halos in 
which they reside, as expected in the CDM paradigm, 
then this implies that even at fixed galaxy mass higher 
a galaxies live in more clustered halos than low a galax- 
ies. The main determinant of a halo's clustering ampli- 
tude is its mass, but it has also been shown in simulations 
that halos with different formation histories also cluster 
differently, with older m ore centrally concentrated halos 
clust e ring more strongly dGao et aU2 005; Wechslcr et al. 
200 1 iWetzel et al.ll2007MGao fc White! 120071; Uing et al. 
20071 : iLi et alJ 120081 ). This means that galaxies with 
higher a at fixed mass could be preferentially be living in 
either more massive or older more centrally concentrated 
halos. 

When considering how galaxies occupy halos it is also 
important to make the distinction between the galaxy 
that has been formed at the center of a given DM halo, 
the central galaxy, and those that initially formed at the 
centers of other halos and have later been accreted into 
this halo, satellite galaxies. The properties of a central 
galaxy will be largely determined by the formation his- 
tory and properties of its host halo, whereas any effect 
the current host halo has on satellite galaxies is limited 
to the time after it has been accreted. 

In the following sections we address whether our re- 
sult can be explained by a relationship between a and 
halo mass for central galaxies (Section 16. ip or for satel- 
lite galaxies (Section I6.2|) . or by a relationship between 



a and halo formation history (Section 16. 3p 

6.1. Central galaxy halo mass correlations 

Since DM halo mass has the largest effect on cluster- 
ing amplitude, the simplest explanation for our results is 
that a is better correlated with its host halo mass than 
stellar or dynamical mass are. For galaxies living at the 
centers of halos, central galaxies, there is a str ong correla- 
tion between stellar mas s and halo mass fe.g. lYang et al.1 
I2009t IWake et all 120111) , but with a significant scatter 
(jYang et al.1 12009 . l2011h . If the scatter between a and 
halo mass for centrals is substantially lower than the scat- 
ter with M s tar this could result in clustering observations 
much as we have observed. 

It is striking in Figures [3] and 2] that the clustering ra- 
tio between high and low a samples becomes even larger 
as the scale decreases. The clustering amplitude on such 
scales is determined by the separations of galaxies within 
individual halos, so between central and satellite galax- 
ies. Thus these smaller scale clustering measurements 
can provide an additional constraint as to whether a 
tighter correlation between a and halo mass (Mhaio) for 
central galaxies could produce the observed clustering 
results. 

We show in Figure [T0l an illustrative attempt at mod- 
eling the effect of a smaller scatter between a and 
Mha io than between M sta r and M na io using the halo model 
(see I Wake et al.|[201lL for details of our model). We have 
fitted halo occupation distributions (HOD), the mean 
number of central and satellite galaxies as a function of 
halo mass, to the clustering measurements for a series 
of samples of SPS S galaxies limited in stellar mass (i.e. 
M sta r>Mi ow ) fsee IWake et all 120121 for details of these 
measurements). Subtracting two of the resulting HODs 
with slightly different Mi ow gives the HOD for a sample of 
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Figure 8. The dependence of large scale bias on elliptical probability (top) and g — r color (bottom) at fixed Md yn (left)and 
a (right). The bias does depended on elliptical probability at fixed mass but not at fixed a. However, there is a significant 
dependence of the bias on color at both fixed mass and a. 



galaxies with a narrow range in M sta r- We show an HOD 
for such a sample as the black line in the top panel of 
Figure [TUl which has been chosen as it matches the low- 
est M s tar sample we use (11.04 < log(M sta r) < 11.19). 
We can then attempt to mimic the effect of splitting into 
high and low a samples under the assumption that a for 
central galaxies is better correlated with halo mass than 
Mgtar is, so that the high a sample will have higher av- 
erage halo masses than the low a sample. We do this by 
slightly adjusting the central galaxy halo mass thresh- 
olds that define the high and low halo mass cut offs in 
our HODs. The HODs of these two samples, representing 
high and low a quartiles at fixed M sta r > are shown as the 
red and blue lines in the top panel of Figure [10] These 
mock high and low a samples each contain 25% of the 
central and satellite galaxies of the full sample, with the 
size of the central HOD controlled by the mass thresh- 
old adjustments and the satellite just random sampled 
to 25% of the number density, i.e. the form of the satel- 
lite HOD is left unchanged. The resulting u> p (r p )ratio 
is shown in the bottom panel of Figure [10] for both the 
auto- and cross-correlations. The cross-correlation is cal- 
culated with an HOD suitable for the parent sample we 
use throughout this work (i.e. log(M sta r) > 10.777). 

The model w p (r p ) ratio shown in Figure [TU] is quite 
similar to the measured ratios shown in Figures [3] [4] and 
[5] showing a constant increase in the cross-correlation 
amplitude on intermediate and large scales and then an 



increasing amplitude at small scales. It is always tempt- 
ing to interpret changes in small scale (1-halo) cluster- 
ing as reflecting changes in the satellite population, but 
this is a clear example of how this is not necessarily 
the case. The satellite population has not changed be- 
tween the mock low and high a samples, and it is the 
change in the central galaxies which is causing a change 
in the numbers of central-satellite pairs that causes the 
increased clustering amplitude on small scales. More ex- 
plicitly, at a given stellar mass the higher a galaxies are in 
more massive halos, which typically contain more satel- 
lite galaxies. This results in an increase in the number 
of central-satellite pairs, thus increasing the small scale 
clustering amplitude. The converse is true for lower dis- 
persion centrals i.e. they are in less massive halos, with 
fewer satellites and so fewer central-satellite pairs. The 
satellite-satellite pairs remain unchanged, but since for 
these massive galaxies the satellite fraction is low so there 
may only be at most one or two satellites in a given halo 
the central-satellite pairs make up a significant contribu- 
tion to the small scale clustering amplitude. 

This exercise clearly demonstrates that if there is a 
much tighter correlation between a and halo mass than 
between M s t ar and halo mass then the clustering ratios 
we measure would be expected. There is, of course, some 
theoretical expectation that this would be the case. The 
stellar mass of a central galaxy depends on a large num- 
ber of complicated astrophysical baryonic processes such 
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Figure 9. The dependence of large scale bias on a at fixed Md yn (left) and Md yn at fixed a (right) for elliptical galaxies (P e i > 
0.6; top) and red galaxies (g — r> 0.9; bottom). The same result as for the whole sample remains, there is significant dependence 
of the bias on a at fixed Md yn but none on Md yn at fixed a. 



as gas cooling and feedback as well as mergers from satel- 
lite galaxies, which may give rise to a large scatter be- 
tween halo mass and stellar mass, cr, on the other hand, 
is tracing the depth of the potential at the center of the 
halo and whilst this may be modified by the evolution of 
the baryons it is likely to be much more directly linked 
to the halo properties and in particular its mass. 

6.2. Satellite galaxy halo mass correlations 

Whilst our example halo model shows that our results 
can be produced just by central galaxies, it is of course 
possible to produce similar clustering trends by modi- 
fying the satellite distribution or both the central and 
satellite distributions. To increase the clustering ampli- 
tude of galaxies at fixed mass by modifying the satellites 
requires that higher a satellites are preferentially found 
in more massive halos. 

One possible process that could create such an effect 
is the stripping of the outer regions of a satellite galaxy 
by tidal interactions as it orbits it parent halo. This 
would likely have a minimal effect on a whilst reducing 
the mass and size of the galaxy. The likelihood of a 
galaxy getting disrupted in this way increases as the halo 
mass increases but decreases as the mass of the galaxy 
increases. Since, the galaxies in this study are all more 
massive than 6 x 1O 1O M0 it would only be effective for 
satellites passing very close to the center of the most 
massive halos an occurrence that would be quite rare. 



However, assuming that this process was an effective 
one we would expect to see its affects on the cluster- 
ing both when the mass is held fixed and a varied and 
visa-versa. Satellite galaxies at fixed mass would have 
higher a in higher mass halos, since they would have 
been stripped more. This would mean that the satellite 
galaxies in the higher a samples would be preferentially 
found in higher mass halos and hence would cluster more 
strongly, just as we have observed. The opposite would 
be true of satellite galaxies with fixed tr; one would expect 
them to have lower masses in more massive halos since 
again they would have experienced more tidal stripping. 

One could then imagine a scenario that could explain 
our observations. If we assume that the more massive 
a central galaxy or the higher its a then the more mas- 
sive its parent halo, then for centrals one would observe 
higher clustering amplitudes both for higher mass galax- 
ies at fixed a and higher a galaxies at fixed mass. At 
fixed a the higher mass satellites would typically be in 
lower mass halos (where they haven't been stripped) and 
visa-versa, resulting in an overall reduction in the clus- 
tering ratio between the high and low mass samples. The 
opposite would be true at fixed mass, where the higher 
a satellite galaxies would preferentially be in higher mass 
halos (where the stripping is effective), increasing the 
clustering ratio between high and low a samples. A sim- 
ilar effect could also occur when considering Md yn rather 
than M st ar since the size of a galaxy would be reduced 
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Figure 10. Top: An example halo occupation distribution 
representing the case where correlation between a and dark 
matter halo mass is tighter than the correlation between 
M star and halo mass. The black line represents an HOD for 
a galaxy population with a narrow slice in M sta r. The red 
and blue lines each contain 25% of the M sta rslice sample and 
have had their central halo mass thresholds adjusted to sam- 
ple slightly higher and lower halo masses from within that 
sample, to represent high and low a samples. In each case 
the dashed lines show the central galaxies and the dotted the 
satellites. Bottom: The auto- (solid) and cross- (dotted) cor- 
relation function ratios between the example high and low 
aHODs. They show the very similar characteristics to the 
measured correlation function ratios on both small and large 
scales. 

as it is stripped, thus reducing Md yn - This does seem to 
be a viable explanation of our results, although we note 
that the fraction of satellite galaxies should be relatively 
low in our samples, < 20%. 

6.3. Central galaxy halo age correlations 

Finally, as already mentioned, there is an alternative 
to the halo mass dependence outlined above, whereby in- 
stead of a tight correlation between central galaxy a and 
halo mass there is a residual correlation between a and 
a halo's formation history or age. N-body simulations 
show that at a given mass halos that formed earlier 
(and are more concentrated) ar e more clustered than 
younger hal os (iGao et al.l 120051: IWechsler et al. 2006; 
Wetzel et all [20071 : iGao fc White! l2007t Uing et aill2007t 
Li et al.ll2008h . This is an appealing explanation since it 
would link together both the correlation between a and 
clustering amplitude as well as that between a and the 
star formation history of galaxies. At a given mass galax- 
ies with higher a exhibit older stellar populations which 



Figure 11. Top: The relationship between R e and a at hxed 
M s tar (10.85 < log(M atar ) < 10.9). The blue points show the 
running median. At a given M star galaxies with higher a are 
smaller. Bottom: The relationship between M star and a at 
fixed M d yn (11.35 < log(M dyn ) < 11.45), showing very little 
trend. 

would be naturally explained by the earlier formation of 
their parent halos and thus their higher clustering am- 
plitude. 

6.4. Consequences of these relationships 

Assuming that a tighter relationship exists between 
a and halo properties than between galaxy mass and 
halo properties for central galaxies, as outlined in Sec- 
tions 16.11 and 16.31 we can use our measurements to infer 
how a, M s tar, Md yn and R e interrelate with one another, 
and with halo mass and concentration. We expect that 
it, M s tar and Mdyn all depend on halo mass and each 
other as follows: 



GM. 



dyn 



halo 



Re 



R e 



Rvir /c 



(5) 



where M na io is the halo mass, R V i r i s the halo virial radiu s 
and c is the halo concentration ((Hopkins "etall 12009ft . 
Therefore, when M sta r is fixed and a varied either R e , 
Mhaio or c must vary. Since the clustering amplitude does 
change then either M n aio or c must also change, such that 
higher dispersion galaxies must live in more massive or 
more concentrated (older) halos. Whilst it is the case 
that galaxies with higher a at a fixed M sta r typically 
have a smaller R e (see Figure [11]) they must also occupy 
more massive or more concentrated halos. Conversely 
as M s tar is increased at fixed a R e must be changing 
as well such that more massive galaxies have larger R ei 
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since there is no change in the clustering amplitude so 
Mhaio and c aren't varying. This means that at fixed 
a increasing M sta r increases the ratio between M sta r and 

M ha l . 

We see a similar effect with Md yn - When Md yn is held 
fixed and a varied R e must also vary such that higher 
a leads to lower R e . We also know that either Mhaio or c 
increases with a at fixed Md yn since the clustering ampli- 
tude increases. M sta r could potentially vary in either di- 
rection with the size of the change depending on the ratio 
of the contributions of dark and stellar mass to the poten- 
tial and the magnitude of the change in Mhaio or c with a. 
In fact Figure [TT] shows that at fixed Md yn there is very 
little average change in M star as a increases, indicating 
that Mhaio or c have to increase with a at fixed Md yn just 
as we observe. When Md yn is varied at fixed a R e must 
also change, increasing as Md yn increases. Since the clus- 
tering amplitude does not change, Mhaio and c do not 
change and so M sta r must be increasing. This implies 
that at fixed a galaxies with higher Md yn will have both 
larger M star and R e . 

7. SUMMARY AND CONCLUSIONS 

We have investigated how the clustering of galaxies de- 
pends on stellar mass, dynamical mass, velocity disper- 
sion, surface mass density, color and morphology using 
large samples of massive galaxies from the SDSS. We 
split the samples into narrow ranges in mass, a or E and 
measure any residual clustering dependence on the other 
parameters. We find the following: 

1. When M star or Md yn are fixed there is a signif- 
icant dependence of the clustering amplitude on 
a on all scales, such that galaxies with higher a are 
more strongly clustered. Conversely when a is fixed 
there is no dependence of the clustering amplitude 
on either M sta r or Md yn - These trends remain when 
we limit the samples to galaxies that are red or 
morphologically elliptical. We see a similar trend 
using S instead of a. 

2. When M sta r or Md yn are fixed there is a weak de- 
pendence of the clustering amplitude on morphol- 
ogy. When a is fixed there is no significant mor- 
phological dependent clustering. 

3. There always remains a strong dependence of the 
clustering amplitude on a galaxy's g ~ r color at 
fixed mass or a. This is most likely caused by satel- 
lite galaxies in massive halos always being prefer- 
entially redder than central galaxies of the same 
mass or a. 

We suggest that for our main finding, point 1. above, 
that there are three possible explanations; the relation- 
ship between a or E and halo mass for central galax- 
ies is tighter than that between M sta r or Md yn and halo 
mass; the relationship between a or E and halo age (or 
concentration) for central galaxies is tighter than that 
between M sta r or Md yn and halo age (or concentration); 
tidal stripping of satellite galaxies reduces the size and 
mass of the galaxy whilst having a minimal effect on 
a and is more effective in more massive halos. Or it 
could be a combination of all three effects. If the host 
halo mass or age are indeed better correlated with a and 



E than with M st ar, these halo properties may also drive 
the star formation hist ory of galaxies, explaining why 
Kauffman n et al.l (|2003T ) found that star for mation rate is 
better correlated with E than M sta r, and why lFranx et al.1 
(2008) found that color and star formation rate are better 
correlated with a than M sta r ■ 
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APPENDIX 
A. SAMPLE PROPERTIES 



Table 3 

Details of the high and low stellar mass samples at fixed velocity dispersion 



Velocity Dispersion (km/s) log(Stcllar Mass) (Mq) 

Min Max Mean Median Min Max Mean Median Ngal 



210.0 


242.3 


223.5 


222.5 


10.78 


10.91 


10.85 


10.85 


3242 


210.0 


242.3 


226.4 


226.3 


11.16 


11.93 


11.30 


11.26 


3242 


242.3 


274.6 


255.2 


253.8 


10.78 


10.98 


10.89 


10.89 


1619 


242.3 


274.6 


257.8 


257.5 


11.27 


11.88 


11.40 


11.37 


1619 


274.5 


306.4 


286.3 


284.5 


10.78 


11.07 


10.95 


10.96 


591 


274.5 


306.8 


288.6 


287.2 


11.37 


11.91 


11.50 


11.46 


590 



Table 4 

Details of the high and low velocity dispersion samples at fixed stellar mass 



log(Stcllar Mass) (Mq) Velocity Dispersion (km/s) 

Min Max Mean Median Min Max Mean Median Ngal 



11.04 


11.19 


11.10 


11.09 


79.0 


177.2 


155.5 


159.3 


3242 


11.04 


11.19 


11.12 


11.12 


232.5 


921.7 


259.7 


252.3 


3242 


11.19 


11.33 


11.24 


11.24 


92.6 


198.4 


173.6 


178.0 


1601 


11.19 


11.33 


11.26 


11.26 


253.6 


565.1 


278.5 


273.0 


1600 


11.33 


11.48 


11.39 


11.38 


108.3 


219.6 


192.3 


198.3 


610 


11.33 


11.48 


11.40 


11.40 


273.7 


562.7 


297.3 


291.8 


610 



Table 5 

Details of the high and low dynamical mass samples at fixed velocity dispersion 
Velocity Dispersion (km/s) log(Dynamical Mass) (Mq) ' 



Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


210.0 


242.3 


222.4 


221.0 


10.67 


11.21 


11.12 


11.13 


3242 


210.0 


242.3 


227.0 


227.2 


11.45 


12.22 


11.59 


11.55 


3242 


242.3 


274.6 


254.7 


253.0 


10.51 


11.33 


11.23 


11.24 


1619 


242.3 


274.6 


258.3 


257.7 


11.60 


12.20 


11.73 


11.70 


1619 


274.6 


306.4 


286.1 


284.0 


10.66 


11.44 


11.33 


11.34 


591 


274.6 


306.8 


289.0 


287.9 


11.72 


12.28 


11.86 


11.82 


590 



Table 6 

Details of the high and low velocity dispersion samples at fixed dynamical mass 



log(Dynamical Mass) (Mq) Velocity Dispersion (km/s) 

Min Max Mean Median Min Max Mean Median Ngal 



11.31 


11.47 


11.37 


11.36 


81.4 


183.2 


162.2 


166.3 


3242 


11.31 


11.47 


11.39 


11.39 


235.0 


454.5 


258.5 


252.7 


3242 


11.47 


11.62 


11.53 


11.52 


102.9 


205.9 


183.0 


187.6 


1647 


11.47 


11.62 


11.55 


11.54 


257.4 


410.7 


280.7 


275.4 


1646 


11.62 


11.78 


11.68 


11.68 


115.9 


231.9 


207.0 


213.1 


695 


11.62 


11.78 


11.70 


11.70 


276.7 


500.7 


298.2 


292.7 


695 
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Table 7 

Details of the high and low stellar mass samples at fixed surface stellar mass density 



log(Surfacc mass density) (M@Kpc ^) log(Stcllar Mass) (Mq) 



Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


9.00 


9.16 


9.08 


9.08 


10.78 


10.87 


10.82 


10.82 


3541 


9.00 


9.16 


9.08 


9.07 


11.10 


11.94 


11.24 


11.21 


3541 


9.16 


9.32 


9.24 


9.24 


10.78 


10.86 


10.82 


10.82 


2272 


9.16 


9.32 


9.24 


9.23 


11.08 


12.07 


11.22 


11.18 


2272 


9.32 


9.48 


9.39 


9.39 


10.78 


10.84 


10.81 


10.81 


1075 


9.32 


9.48 


9.39 


9.38 


11.05 


11.75 


11.18 


11.15 


1075 


9.48 


9.64 


9.55 


9.54 


10.78 


10.84 


10.80 


10.81 


383 


9.48 


9.64 


9.55 


9.54 


11.04 


11.78 


11.17 


11.14 


383 



Table 8 

Details of the high and low surface stellar mass density samples at fixed stellar mass 



log(Stcllar Mass) (Mq) 


log(S 


urfacc mass density) (MqK'pc A ) 




Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


10.80 


10.87 


10.83 


10.83 


6.58 


8.71 


8.46 


8.53 


3541 


10.80 


10.87 


10.83 


10.83 


9.15 


10.51 


9.33 


9.29 


3541 


10.87 


10.93 


10.90 


10.90 


6.58 


8.73 


8.48 


8.56 


2761 


10.87 


10.93 


10.90 


10.90 


9.15 


10.24 


9.32 


9.29 


2761 


10.93 


10.99 


10.96 


10.96 


6.77 


8.74 


8.49 


8.57 


2377 


10.93 


10.99 


10.96 


10.96 


9.14 


10.41 


9.31 


9.28 


2377 


10.99 


11.05 


11.02 


11.02 


7.05 


8.74 


8.51 


8.58 


1963 


10.99 


11.05 


11.02 


11.02 


9.14 


10.37 


9.31 


9.27 


1963 



Table 9 

Details of the high and low elliptical probability samples at fixed dynamical mass 



log(Dynamical mass) (Mq) Elliptical probability 



Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


11.31 


11.47 


11.38 


11.37 


0.00 


0.29 


0.09 


0.07 


3911 


11.31 


11.47 


11.39 


11.39 


0.84 


1.00 


0.94 


0.94 


2011 


11.47 


11.62 


11.53 


11.52 


0.00 


0.29 


0.09 


0.06 


1405 


11.47 


11.62 


11.54 


11.54 


0.84 


1.00 


0.94 


0.94 


1869 


11.62 


11.78 


11.69 


11.68 


0.00 


0.29 


0.09 


0.07 


364 


11.62 


11.78 


11.70 


11.69 


0.84 


1.00 


0.95 


0.95 


1272 



Table 10 

Details of the high and low elliptical probability samples at fixed velocity dispersion 



Velocity dispersion (km/s) Elliptical probability 



Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


210.0 


242.3 


223.1 


221.7 


0.00 


0.76 


0.35 


0.37 


4026 


210.1 


242.3 


226.1 


225.9 


0.86 


1.00 


0.95 


0.94 


2471 


242.3 


274.6 


255.3 


254.1 


0.00 


0.77 


0.40 


0.43 


1169 


242.3 


274.7 


256.4 


255.4 


0.86 


1.00 


0.95 


0.95 


1925 


274.7 


306.4 


286.6 


284.9 


0.02 


0.74 


0.44 


0.47 


307 


274.7 


306.8 


287.6 


286.0 


0.86 


1.00 


0.95 


0.95 


871 



Table 11 

Details of the high and low g-r color samples at fixed dynamical mass 



log(Dynamical mass) 


(Mq) 




g- 


r color 






Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


11.31 


11.47 


11.38 


11.37 


0.00 


0.91 


0.75 


0.76 


3260 


11.31 


11.47 


11.39 


11.38 


0.96 


1.87 


0.99 


0.98 


3203 


11.47 


11.62 


11.53 


11.52 


0.00 


0.95 


0.79 


0.81 


1628 


11.47 


11.62 


11.54 


11.53 


0.97 


1.79 


1.00 


0.99 


1660 


11.62 


11.78 


11.69 


11.68 


0.53 


0.96 


0.85 


0.87 


713 


11.62 


11.78 


11.70 


11.69 


0.98 


1.95 


1.01 


1.00 


680 
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Table 12 

Details of the high and low g-r color samples at fixed velocity dispersion 



Velocity dispersion (km/s) 




g 


-r color 






Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


210.0 


242.3 


222.8 


221.3 


0.00 


0.94 


0.87 


0.89 


3280 


210.0 


242.3 


225.9 


225.9 


0.97 


1.95 


1.00 


0.99 


3228 


242.4 


274.6 


254.7 


252.9 


0.00 


0.95 


0.91 


0.92 


1643 


242.3 


274.6 


257.1 


256.6 


0.98 


1.49 


1.01 


1.00 


1606 


274.6 


306.6 


286.1 


284.5 


0.44 


0.96 


0.93 


0.94 


583 


274.7 


306.9 


288.0 


286.6 


0.99 


1.50 


1.02 


1.01 


593 



Table 13 

Details of the high and low dynamical mass samples at fixed velocity dispersion and elliptical probability > 0.6 



Velocity Dispersion (km/s) log(Dynamical Mass) (Mq) 



Min 


Max 


Mean 


Median 


Min 


Max 


Mean 


Median 


Ngal 


210.3 


230.8 


219.3 


218.8 


10.68 


11.20 


11.11 


11.12 


1483 


210.3 


230.8 


221.4 


221.7 


11.43 


12.16 


11.56 


11.52 


1483 


230.8 


251.3 


239.4 


238.7 


10.51 


11.28 


11.19 


11.20 


1160 


230.9 


251.3 


241.9 


242.1 


11.53 


12.19 


11.67 


11.63 


1160 


251.4 


271.8 


260.0 


259.3 


10.88 


11.35 


11.26 


11.26 


765 


251.4 


271.8 


261.3 


260.9 


11.62 


12.18 


11.75 


11.71 


765 



Table 14 

Details of the high and low velocity dispersion samples at fixed dynamical mass and elliptical probability > 0.6. 



log(Dynamical Mass) (Mq) Velocity Dispersion (km/s) 

Min Max Mean Median Min Max Mean Median Ngal 



11.31 


11.45 


11.37 


11.36 


132.0 


206.8 


189.7 


193.1 


1483 


11.31 


11.45 


11.39 


11.39 


246.7 


454.5 


268.3 


261.9 


1483 


11.45 


11.59 


11.51 


11.50 


136.0 


221.0 


204.2 


207.8 


1030 


11.45 


11.59 


11.52 


11.52 


261.4 


386.6 


283.9 


278.5 


1030 


11.59 


11.73 


11.65 


11.64 


156.9 


237.6 


220.6 


223.7 


574 


11.59 


11.73 


11.66 


11.66 


277.2 


500.7 


298.0 


292.2 


574 



Table 15 

Details of the high and low dynamical mass samples at fixed velocity dispersion and g — r > 0.9 
Velocity Dispersion (km/s) log(Dynamical Mass) (Mq) 



Min Max Mean Median Min Max Mean Median Ngal 



210.0 


231.5 


219.1 


218.4 


10.73 


11.19 


11.11 


11.12 


1973 


210.0 


231.5 


221.6 


221.9 


11.42 


12.16 


11.55 


11.51 


1972 


231.5 


253.0 


240.3 


239.5 


10.71 


11.27 


11.18 


11.19 


1442 


231.5 


252.9 


243.0 


243.3 


11.53 


12.19 


11.66 


11.63 


1442 


253.0 


274.4 


261.9 


261.3 


10.96 


11.34 


11.25 


11.25 


872 


253.0 


274.4 


263.6 


263.7 


11.62 


12.20 


11.75 


11.71 


872 



Table 16 

Details of the high and low velocity dispersion samples at fixed dynamical mass and g — r > 0.9. 
log(Dynamical Mass) (Mq) Velocity Dispersion (km/s) 



Min Max Mean Median Min Max Mean Median Ngal 



11.31 


11.46 


11.37 


11.36 


125.5 


203.8 


187.0 


190.3 


1973 


11.31 


11.46 


11.39 


11.39 


245.4 


454.5 


267.0 


260.7 


1972 


11.46 


11.60 


11.52 


11.51 


129.3 


220.0 


202.4 


206.1 


1213 


11.46 


11.60 


11.53 


11.53 


262.3 


410.7 


284.9 


279.0 


1212 


11.60 


11.75 


11.66 


11.65 


159.5 


237.7 


219.9 


223.5 


624 


11.60 


11.75 


11.67 


11.67 


278.4 


500.7 


299.5 


293.8 


624 
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Figure 12. The redshift, velocity dispersion, stellar mass and elliptical probability distributions for the high (red) and low 
(blue) stellar mass at fixed velocity dispersion samples (left) and the high and low velocity dispersion at fixed stellar mass 
samples (right) The solid and dashed vertical lines in the velocity dispersion and stellar mass plots show the mean and median 
of the distributions respectively. 
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Figure 13. The redshift, velocity dispersion, dynamical mass and elliptical probability distributions for the high (red) and low 
(blue) dynamical mass at fixed velocity dispersion samples (left) and the high and low velocity dispersion at fixed dynamical 
mass samples (right). The solid and dashed vertical lines in the velocity dispersion and stellar mass plots show the mean and 
median of the distributions respectively. 
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Figure 14. The redshift, surface mass density, stellar mass and elliptical probability distributions for the high (red) and low 
(blue) stellar mass at fixed surface mass density samples (left) and the high and low surface mass density at fixed stellar mass 
samples (right). The solid and dashed vertical lines in the surface mass density and stellar mass plots show the mean and 
median of the distributions respectively. 
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Figure 15. The redshift, g — r color, dynamical mass and elliptical probability distributions for the high (red) and low (blue) 
g — r color at fixed dynamical mass samples (left). The redshift, velocity dispersion, g — r color, and elliptical probability 
distributions for the high (red) and low (blue) g — r color at fixed velocity dispersion samples (right). The solid and dashed 
vertical lines show the mean and median of the distributions respectively. 
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Figure 16. The redshift, velocity dispersion, dynamical mass and elliptical probability distributions for the high (red) and low 
(blue) elliptical probability at fixed dynamical mass density samples (left) and at fixed velocity dispersion samples (right). The 
solid and dashed vertical lines show the mean and median of the distributions respectively. The distributions in the redshift 
and dynamical mass plots have been scaled to the same number of objects. 
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B. w p (r p ) RATIOS FOR SAMPLES CUT BY ELLIPTICAL PROBABILITY OR COLOR 
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Figure 17. The ratio of the projected correlation functions between high and elliptical probability samples at fixed dynamical 
mass (left) and velocity dispersion (right). The best fit ratio on scales 1.6 < r p < 25.1 h _1 Mpc is shown as the dashed line. 
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Figure 18. The ratio of the projected correlation functions between red and blue g-r color samples at fixed dynamical mass 
(left) and velocity dispersion (right). The best fit ratio on scales 1.6 < r p < 25.1 h~ 1 Mpc is shown as the dashed line. 
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Figure 19. The ratio of the projected correlation functions between high and low a samples at fixed Md yn for Elliptical galaxies 
(P e i > 0.6 left) and between high and low Md yn samples at fixed a for Elliptical galaxies (P e ; > 0.6 left). The best fit ratio on 
scales 1.6 < r p < 25.1 h _1 Mpc is shown as the dashed line. 
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Figure 20. The ratio of the projected correlation functions between high and low a samples at fixed Md yn for red galaxies 
(g — r > 0.9 left) and between high and low Md yn samples at fixed a for red galaxies (g — r > 0.9 left). The best fit ratio on 
scales 1.6 < r p < 25.1 h _1 Mpc is shown as the dashed line. 
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